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^ Abstract 

S The two-dimensional multifractal detrended fluctuation analysis is applied to reveal the multifractal properties of the fracture surfaces of foamed 
, ; polypropylene/polyethylene blends at different temperatures. Nice power-law scaling relationship between the detrended fluctuation function 
d Fq and the scale s is observed for different orders q and the scaling exponent h[q) is found to be a nonlinear function of q, confirming the 
Cj ()resence of multifractality in the fracture surfaces. The multifractal spectra /(a) are obtained numerically through Legendre transform. The 
shape of the multifractal spectrum of singularities can be well captured by the width of spectrum Aa and the difference of dimension A/. 
With the increase of the PE content, the fracture surface becomes more irregular and complex, as is manifested by the facts that Aq increases 
and A/ decreases from positive to negative. A qualitative interpretation is provided based on the foaming process. 
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1. Introduction 

' The plastic foam industry is fast growing and the plastic 
foams have drawn a great deal of interest in recent decades yj]. 
The cell sizes of the plastic foams strongly influence the ap- 
plication of the materials IZO. The microcellular polymer, with 
cell sizes hundreds of times smaller than those of conventional 
plastic foams, can offer some unique properties that conven- 
tional foams do not possess, such as higher impact strength, 
higher toughness, higher stiffness-to-weight ratio, higher fa- 



tigue life, higher thermal stability, lower dielectric constant, 
and lower thermal conductivity yl Ij, jsl |6|] . Therefore, micro- 
cellular polymer has wide industrial and everyday-life applica- 
tions including food packaging, airplane and automotive parts, 
sporting equipment, insulation, controlled release devices and 
filters, and so on ITl |8|]. The sizes and morphology of cells 
have considerable influence on the application of foaming ma- 
terials, which are usually studied using the fracture surfaces 
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lauiiH]. 

In most cases, fracture surfaces of different materials are 
self-similar, which can be characterized by fractal and multi- 
fractal theories tlSJ. The monofractal properties of the mor- 
phology of fracture surfaces have been investigated for metals 
liSInLceramics ifliElH, poWmersjEiySElJii 
cretes 25U2S|23|2i, alloys IM S laJIiiWlsr 
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371 bsl 13911 . and many other materials. Moreover, the 
multifractal features of surfaces have also been studies. Raoufi 
et al. analyzed the multifractal spectrum of ITO thin films pre- 
pared by electron beam deposition method and the /(a) shapes 
of ITO thin films remained left hooked after annealing at 200° C 
and 300°C iH. Moktadir et al. researched the multifractal 
properties of Pyrex and silicon surfaces blasted with sharp par- 
ticles, and found that the long-range correlations were the ori- 
gin of the multifractal behaviour 114 ill . 

In this work, we investigate the multifractal properties of 
the fracture surfaces of polypropylene (PP) and polyethylene 
(PE) blends foamed with supercritical carbon dioxide. The two- 
dimensional multifractal detrended fluctuation analysis (MF- 
DFA) is adopted, which has the advantages of easy implemen- 
tation, high precision, and low computational time [Al]. The 
MF-DFA approach has been applied to investigate the land- 
scape of the Yardangs region on Mars and the fracture sur- 
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Nomenclature 


s 


scale of boxes 


X{i,j) 


two-dimensional matrix 


Uv.wiiJ) 


cumulative sum 


Uv,w{i,j) 


fitting bivariate polynomial 


^v,w{hj) 


residual matrix 


Q 


order of detrended fluctuation function 


F,{s) 


detrended fluctuation function 


Df 


fractal dimension 


h{q) 


scaling exponent function 


r{q) 


mass exponent function 


a{q) 


singularity strength function 


Q^max 


maximum singularity 


Q^min 


minimum singularity 


Aa 


width of multifractal spectrum 


/(«) 


multifractal singularity spectrum 


A/ 


difference. A/ = /(ofmax) - /(amin) 


T 


foaming temperature 


w 


weight fraction of PE 


P 


probabiUty measure 



face of a foamed polyurethane sample with supercritical carbon 
dioxide 114211 . the combustion flames in four-burner impinging 
entrained-flow gasifier ll43ll and Pollocks's drip paintings [441. 



2. Experimental 

2.1. Materials 

The isotactic polypropylene (Y1600) we have used was pur- 
chased from Shanghai Petrochemical Company, China. The 
crystallinity and melting temperature of the isotactic polypropy- 
lene was 47% and 169°C, respectively. 

The low density polyethylene (2426H) was purchased from 
Yangzi Petrochemical Company, China. The crystallinity and 
melting temperature of the low density polyethylene was 42% 
and 111.2°C, respectively. 

The CO2 (purity: 99.9%) supplied by Shanghai Air Product 
Company, China, was utilized as a blowing agent. 



an integrated backflow channel. The blending was carried out 
for 10 minutes under a 0.6 MPa nitrogen atmosphere with a 
blending temperature of 190° C and a screw speed of 50 rounds 
per minute. After blending, the rod-like PP/PE blends were col- 
lected at the die exit for foaming. 

Fig. [U illustrates the schematic experimental setup of the 
depressurization batch foaming process. A high-pressure vessel 
in stainless steel was used. The internal volume of the vessel 
was 80 cm'^, calibrated with distilled water by a syringe pump. 
A pressure transducer of type P31 from Beijing Endress & 
Hauser Ripenss Instrumentation Company Limited, was used 
to measure the pressure with a precision of ±0.01 MPa and a 
valve of type Swagelok SS-1RS8MM to release the CO2 gas. 
A computer installed with a PCI bus data acquisition system 
was connected to the above pressure transducer to record the 
pressure decay during a depressurization process. 
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Fig. 1. Schematic diagram of the experiment setup for the batch foaming 
process: 1, CO2 steel cylinder; 2, syringe pump; 3, oil bath; 4, high-pressure 
vessel; 5, stirrer; 6, temperature controller; 7, pressure transducer; 8, data 
acquisition system. 

In the foaming process, all PP/PE blends samples (including 
the pure PP sample for comparison) were placed in the high- 
pressure vessel to ensure the same foaming condition and the 
vessel was purged with low-pressure CO2. Thereafter, a given 
amount of CO2 was charged. The CO2 loading was achieved by 
a DZB-IA syringe pump of Beijing Satellite Instrument Com- 
pany (China) with a precision of 0.01 cm'^. The high-pressure 
vessel was immersed in a silicone oil bath and rapidly heated 
to a desired saturation temperature. After the sorption of CO2 
into the blends samples reached a sort of equilibrium, the CO2 
was released rapidly from the high-pressure vessel. The foamed 
blends samples were taken out for subsequent analysis. The 
samples were immersed in liquid nitrogen for 10 minutes and 
then fractured. The cell morphologies of the foamed blends 
samples were characterized by a JSM-6360LV scanning elec- 
tron microscopy (SEM). The polymer foaming was influenced 
by many factors. In this work, we consider the influence of dif- 
ferent mixture ratios at two different temperature T = 150 and 
140°C. 



2.2. Experiment process 



2.3. SEM images 



Five different weight fractions of the PP/PE blends were 
studied, in which the PP fractions are 100%, 95%, 90%, 75% 
and 50%. The blends were prepared in a Haake Minilab sys- 
tem, which is based on a conical twin-screw compounder with 



From the experiments, many SEM pictures of the fracture 
surface of the foams under different experimental conditions 
were obtained. Fig.|2]illustrates a typical image of fracture sur- 
face of the foamed sample prepared at 150°C and 25 MPa with 



a depressurization rate of 200 MPa/s and the mixture ratio PP : 
PE = 90 : 10. The depressurization rate was not a constant 
during the depressurization process. The rate of 200 MPa/s was 
the largest depressurization rate during the whole depressuriza- 
tion process, which is the key parameter mainly determining 
the cell nucleation rate. The data acquisition system obtained 
the change of pressure with time and then we were able to 
calculate the depressurization rate at any time, such that the 
depressurization rate can be calculated as the ratio of the pres- 
sure difference over the depressurization time. The size of the 
images was 960 x 1280 pixels, and we intercepted them into 
800 X 1200 to eliminate the noise of the mark on the images. 
The images were stored in the computer as two-dimensional 
arrays in 256 grey levels for multifractal analysis. 




Fig. 2. SEM micrograph of the foamed sample (PP : PE = 90 ; 10) under 
condition: 150° C, 25 MPa, depressurization rate 200 MPa/s. 



3. Two-dimensional multifractal detrended fluctuation 
analysis 

The multifractal properties may exist in the fracture surface 
of the foams. To unveil the multifractal characteristics of the 
SEM images, the two-dimensional MF-DFA is applied. This 
method has become popular for its simplicity and easy com- 
puter implementation. The two-dimensional MF-DFA is the ex- 
tension of the DFA method which is frequently used in com- 
puting the roughness exponent of monofractal signals and the 
identification of long range correlations in non-stationary time 
series 14511 . The two-dimensional MF-DFA applied to our ex- 
periment data can be summarized as following steps: 

Step 1. Consider a self-similar (or self-affine) surface, which 
is denoted by a two-dimensional array X{i,j), where i = 
1, 2, ...,M and j ~ 1, 2, ..., N. The surface is partitioned into 
Ms X Ns disjoint square segments of the same size sxs, where 
Ms = [^^/s] and Ns = [-^/s]- Each segment can be denoted 
by Xu -u, such that X„..ii,(i,j) ^X{li+i,l2+j)forl ^i,j ^ 
s, where li = {v — l)s and h = {w — l)s. 

Since M and N are often not a multiple of the segment 
size s, two orthogonal strips at the end of the profile may 
remain and some data will be ignored by this way. In order to 
take these ending parts of the surface into consideration, the 
same partitioning procedure can be repeated starting from the 
other three corners. Then there will be AMgNs segments and 
calculating the average over them can eliminate the boundary 
influence. 

Step 2. For each segment Xy^^ identified by v and w, the 
cumulative sum Uy,uu{i,j) is calculated as follows: 



Fig. [3] illustrates a typical fracture surface image of the 
foamed sample prepared at 150°C, 25 MPa, with the depres- 
surization rate of 200 MPa/s, and the mixture ratio PP : PE = 
50 : 50. And we can find the open-cell easily, and the same 
result can be obtained at 140°C. Thus, we can imagine that 
there is a threshold content of PE in the sample at the range 
25% ^ 50% where the cells change to open. 
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Fig. 3. SEM micrograph of the foamed sample (PP : PE = 50 : 50) under 
condition: 150°C, 25 MPa, depressurization rate 200 MPa/s. 



Uv,w{iJ) = X! X! Xv,w{ki,k2) , (1) 

fel=l fe2=l 

where I ^ i,j ^ s. Note that Uv,w itself is a surface. 

Step 3. The trend of the constructed surface Ui,,iu can be 
determined by fitting it with a prechosen bivariate polynomial 
function Uy^^{i,i). The parameters of Uy,^{i,j) can be esti- 
mated easily through the least square method fit to the data in 
each segment Xy^^j. Then, the residual matrix can be obtained. 
The detrended fluctuation function F{v,w,s) of the segment 
Xj, ,„ is defined via the sample variance of the residual matrix: 

^ s s 

F'^{v,w,s) ^ —^^[uy^^,{i,j) -Uy^y,{i,j)]'^ . (2) 

* i=l j=l 

Step 4. The overall detrended fluctuation is calculated by 
averaging over all the segments, and the value of the (/th-order 
fluctuation function is 

** * u=l to=l ) 

where q can take any real value except for q = 0. When q = 0, 
we have 

C M. N. ^ 

I, V—l It? — 1 J 



according to L'Hopital's Rule. 

Step 5. Varying the value of s in the range from s„iin w 6 
to Smax ~ min(Af , N)IA, we can determine the scaling relation 
between the detrended fluctuation function Fq{s) and the size 
scale s, which reads: 



F,{s) 



Mq) 



(5) 



where the exponent h{q) is called the generalized Hurst in- 
dex. The scaling exponent h{q) is a constant for monofractals 
and a nonlinear decreasing function of q for multifractals. For 
positive q values, h{q) describes the scaling behaviour of the 
segments with large fluctuations, whereas for negative q values 
h{q) concerns with small fluctuations. 

In the standard multifractal formalism based on partition 
function 14611 . the multifractal nature is characterized by the 
mass exponents T{q), which is a nonlinear function of q. The 
r(g) function is related to h{q) through 



'(?) = qh{q) - Df 



(6) 



where Df is the fractal dimension of the geometric support of 
the multifractal measure l42l 14511 . According to the Legendre 
transform, we can obtain the singularity strength function a{q) 
and the singularity spectrum /(a) as follows 14611 : 



a = T'{q) = h{q) + qh'{q) 
f{a) ^ qa{q) - r(g) . 



(7) 
(8) 



4. Results 



4. 1 . Multifractal analysis 

As a first step, we perform multifractal analysis on the SEM 
images to check if the fracture surfaces possess multifractal 
nature or not. We take the SEM picture showed in Fig. |2] as 
an example to illustrate the power-law scaling between the de- 
trended fluctuation function Fq{s) and the scale s. Fig.|4]shows 
the fluctuation Fq{s) as a function of s for five different values 
of q in double-logarithmic coordinates. The results for q = —3, 
0, 3 and 6 have been shifted upward by 0.4, 0.8, 1.2 and 1.6 
for clarity. According to the figure, the data points for every 
q fall on a straight line, indicating a perfect power-law scal- 
ing between Fq{s) and s, as expressed in Eq. (|5]l. The scaling 
range is from 6 pixels to 200 pixels for all q values. The fluc- 
tuation functions Fq{s) of other fracture surfaces also exhibit 
nice power-law behaviors. 

According to Eq. Q, the slopes of the straight lines illus- 
trated in Fig. m are the scaling exponent h{q), which can be 
determined by simple linear regressions of In F, (s) against In s 
for different q. Fig. |5] illustrates h{q) as a function of q for 
— 6 ^ g ^ 6 with the same image. It is noteworthy to point out 
that, due to the finite size of the data, the overall fluctuation 
values Fq{s) for too large q are not statistically significant. It is 
evident that h{q) is a decreasing nonlinear function of q. The 
mass exponent function T(q) is also calculated numerically us- 
ing Eq. ^, where Df = 2. The inset of Fig. |5] shows the mass 
exponent function T{q) for — 6 ^ g ^ 6. The nonlinearity of 




Fig. 4. Log-log plots of Fq(s) versus s for five different values of q. The 
lines are the best power-law fits to the data. The plots for q = —3, 0, 3 and 
6 are shifted upward by 0.4, 0.8, 1.2 and 1.6 for clarity. 

h{q) and T{q) confirms that the fracture surface under inves- 
tigation possesses multifractal nature. We note that all other 
fracture surfaces also exhibit multifractaUty. 




Fig. 5. Dependence of the scahng exponent h(q) with respect to the order 
q. The inset shows the mass exponent function T[q). 



4.2. Multifractal spectra 

Theoretically, the geometric support of a multifractal mea- 
sure can be decomposed into interwoven fractal sets, each of 
which is characterized by its singularity strength a. The frac- 
tal dimension of the underlying fractal set associated with a is 
/(a), which is the well-known singularity spectrum or multi- 
fractal spectrum. In this vein, a and f{a) are two most impor- 
tant characteristics in the description of the multifractal. 

We have numerically obtained the values of a{q) and f{a) 
through the Legendre transform from T{q) for all the fracture 
surfaces under investigation. Fig. |6] illustrates the multifractal 
spectra f{a) with respect to a for five samples with different 
ratios for PP and PE at 150°C. The multifractal spectra for dif- 
ferent blends samples exhibit different shapes. With the increase 
of the proportion of polyethylene, the singularity spectrum be- 
comes wider It means that the surface with higher polyethy- 
lene proportion is more irregular, as is clear from the compar- 



b 

^ 



ison between Fig. |2] and Fig. [3] Another intriguing feature in 
Fig. |6]is that the f{a) function for PP = 50% becomes neg- 
ative when a is larger than about 2.6. The negative dimension 
if (a) < 0) was investigated in several experiments such as the 
diffusion-limited aggre gati on 14711 and the energy dissipation 
field of turbulent flows 14811 . The negative dimension describes 
the rarely occurring events ll49ll and one needs an exponentially 
increasing number of samples to observe the subsets with the 
same a value 14811. 
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are the singularity strengthes associated with 
the regions of the sets where the measures are the least and most 
singular, respectively. The corresponding /(amin) and /(amax) 
values reflect the fractal dimensions of the two regions charac- 
terized by a = amin and a = amax- The shape of the multi- 
fractal spectrum f{a) can be captured to a great extent by the 
width of the multifractal spectrum Aa = amax — ctmin and the 
difference of the fractal dimensions A/ = /(c^max) ~ /(amin)- 
We discuss the dependence of Aa and A/ with respect to the 
PP/PE ratio. 

5.1. The dependence of Aa with respect to the PP/PE ratio 

In the formalism of multifractal, amin is related to the maxi- 
mum probability measure by Pmax ~ e"""", where e represents 
the scale approaching zero and it is a small quantity, whereas 
amax is related to the minimum probability measure through 
^min ^ £"">">= . The width Aa can be used to describe the range 
of the probability measures: 



JPn 



(9) 
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a 



The greater the Aa value, the wider the probability distribution 
and the larger the foamed growth probability of the surface. 



Fig. 6. Multifractal spectra of the fracture surfaces of different PP/PE samples 
prepared at T = 150°C. The curves for PP = 75%, 90%, 95%, and 100% 
have been shifted downward by 0.2, 0.4, 0.6 and 0.8 for clarity. 

Fig. |7]plots the multifractal spectra /(a) of the singularity 
a for five samples with different PP/PE ratios at temperature 
T = 140° C. Again, the multifractal spectra for different blends 
samples exhibit different shapes and the singularity spectrum 
becomes wider with the increase of the proportion of polyethy- 
lene. However, no negative dimension is observed in this case. 




Fig. 7. Multifractal spectra of the fracture surfaces of different PP/PE samples 
prepared at T = 140°C. The curves for PP = 75%, 90%, 95%, and 100% 
have been shifted downward by 0.2, 0.4, 0.6 and 0.8 for clarity. 



5. Discussion 

There are several fundamental quantities related to the mul- 
tifractal spectrum. The minimum and maximum singularities 




Fig. 8. The width of spectrum Aa versus the content of PE (w) for two 
temperatures. 

Fig. [8] illustrates the relationship between Aa and the con- 
tent of PE in the sample at two different temperatures. We find 
that the value of Aa increases with augmenting PE accretion. 
This well agrees with the experimental results that the num- 
ber of cells in the sample increases with the increase of the PE 
content, and the fracture surface becomes more irregular We 
also find that the two curves at different temperaments share 
almost the same variation tendency. Since PP is a semicrys- 
talline polymer with a high crystallinity degree, it needs a very 
rigorous condition for PP foaming since gases do not dissolve 
in the crystalline regions 115 Ol pll 15211 . According to Fig. [H 
the Aa value of the pure PP sample is merely 0.089, which 
means that the fracture surfaces of pure PP samples are close 
to monofractal. When PE is added in the sample, the PE with a 
low crystalhzation temperature(Tc) can be molten (the melting 



point of LDPE we used is 111.2°C)in 140°C, and both the cell 
nucleation and bubble growth may take place in the region of 
the molten PE. Therefore, mixing with PE makes the sample 
foaming easier, and we obtain wider multifractal spectra with 
the increase of the PE content while the crystalline regions may 
remain intact. Comparing the two curves in Fig. [8] the curve for 
foaming at T = 150°C is above that of T = 140°C, which is 
consistent with the results that higher temperatures cause larger 
foamed regions with more cells and larger cell size JSSll . 

According to Fig. [S] Aa increases very quickly when the 
content of PE is more than 25%. The cell opening may be the 
reason for this transilient phenomenon. In our experiment, the 
soft sections (PE) form minor and dispersed phases and the hard 
sections (PP) form a major melt matrix. Cell opening can be 
initiated and propagated through well-dispersed soft domains 
that are entrapped between growing adjacent cells. Even though 
these soft domains can become elongated as cells grow (i.e., cell 
walls become thinner), cell opening is most likely to be initiated 
at the weakest cell^wall sections because of the embedded soft 
polymer phases |52]. When the content of PE reaches a certain 
value, the cell may be opening and continuous to form a three- 
dimensional network, which leads to a much wider distribution 
of the singularity of the fracture surface. We can observe easily 
the open cells in Fig. [3] at T = 150°C, and at T = 140°C as 
well. 
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Fig. 9. The difference of fractal dimension A/ versus tlie content of PE {w) 
for two temperatures. 

temperature increase or a pressure decrease; (3) cell growth and 
coalescence, a combination of mass transfer and fluid dynamics. 
In bubble growth and coalescence, the gas (CO2) transferred 
from small bubble (with high pressure) to big bubble (with low 
pressure) to form the foamed structure 1157115811 . which may be 
the origin of long-range correlation. However, in-depth mech- 
anism (mass transfer and fluid dynamics) research of foaming 
is needed to support this surmise. 



5.2. The dependence of Af with respect to the PP/PE ratio 



6. Conclusion 



The parameter A/ is also a very important quantity in the 
multifractal analysis. The /(amax) value reflects the fractal 
dimension of the subset of the minimum growth probability 



with Np 



iV„ 



-f{a„ 



while /(amin) reflects 



that of the maximum probability such that Np^^^^^ = -/Vq^j,, ^ 
g-/(amin) Hence, the A/(a) value can describe the ratio be- 
tween the regions that the probability measure distributes most 
concentrated and most rarified 



Np^^^/Nf 



.A/ 



(10) 



Thus, A/ < means that there are more concentrated regions 
than rarified sites, whereas A/ > means the contrary. 

Fig. m depicts the dependence of A/ with respect to the 
content of PE in the sample at two different temperatures. With 
the increase of the PE content, the A/ value changes from 
positive to negative, which reflects the change of multifractal 
spectrum shape from left-hooked to right-hooked as showed in 
Fig. m and Fig.|2] It implies that the fracture surface becomes 
more complex and singular when there are more PE in the 
sample. 

The observed multifractal behaviour can be partly interpreted 
in terms of the spatial intermittency l54l 15511 and the origin of 
the multifractal hidden in the fracture surface may be the long 
range correlations of the intermittent fluctuation ll4llJ56ll. A ba- 
sic foaming process can be divided into three steps \4\: (1) mix- 
ing, formation of a homogeneous solution composed of foam- 
ing agent and polymer melt; (2) cell nucleation, phase separa- 
tion induced by a thermodynamic instability which is usually a 



Foaming is a very complex process, and it is a challenging 
task to underpin the exact mechanism explaining the nucleation 
and bubble growth 15711 . In this work, we obtained the foamed 
mixed polymers by batch foaming at fixed pressure and de- 
pressurization rate. We used the two-dimensional multifractal 
detrended fluctuation analysis to unveil the multifractal prop- 
erties of the fracture surfaces of formed PP/PE blends samples. 
Perfect power-law scaling is observed and nonlinear relation- 
ship between the scaling exponent h{q) and q in the consider- 
ing moment order —6 ^ q ^ 6 show that the fracture surface 
exhibits multifractal nature. 

The two important multifractal parameters Aa and A/ were 
calculated to describe the multifractal nature of the fracture 
surfaces. We have found that, with the increase of the PE frac- 
tion, the fracture surfaces become more and more irregular and 
complex, which is indicated by the facts that Aa increases and 
A/ decreases. These two parameters can serve as "complexity 
measures" of the fracture surfaces of foamed polymer blends, 
and other surfaces as well. 
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